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. Abstract. Microscopic descriptions of exotic nuclei are the subject of much experimental and theo- 

retical effort. Not only are such important in their own right but are also necessary for applications in 
Q I nuclear astrophysics. Evaluations of model wave functions may be done with analyses of elastic and 

. inelastic scattering from hydrogen. Those require a realistic model of nucleon-nucleus scattering as 

scattering from hydrogen translates to proton scattering in the inverse kinematics. The Melbourne 
^-folding model for intermediate energy is presented along with various examples. Implications for 
(N , existing and future experimental and theoretical work are discussed. 

T— I ■ 

> : 

O INTRODUCTION 

oo . 

O ■ 

O I The microscopic facets of the structures of exotic nuclei have been the subject of in- 

■ creased experimental and theoretical study, particularly with data now available for the 

. scattering of heavy ions from hydrogen at intermediate energies. Such data are comple- 

"^H i mentary to those obtained for breakup reactions, which only probe the asymptotic part 

-'-^ • of the initial state wave function of the exotic nucleus [1,2]. Also, it has been established 

. that the breakup of ^He is a two-step process [3], with the final state interactions greatly 

^ i influencing the reaction process. To investigate the wave functions of exotics at a deeper 

^ ■ level one requires analyses of complementary scattering data. 

>• ■ In the absence of electron scattering data, intermediate-energy proton scattering rep- 

^ ! resents the best probe of the microscopic structures of exotic nuclei. This was illustrated 

I by Karataglidis et al. who compared electron and proton elastic and inelastic scattering 

■ — ' from ^^C [4] and ^^^Li [5]. They found that the behaviour with momentum transfer of 

the form factors from electron scattering was found also in the differential cross sec- 
tions from intermediate energy proton scattering, reflecting the deficiencies in the wave 
functions from the underlying assumed model structure. Such comparisons were only 
possible when credible models of scattering for both electrons and protons from nuclei 
were specified. 

It is well-known that few -body descriptions of exotics, especially the halo nuclei, are 
able to describe breakup reactions as those models are able to give the correct asymptotic 
behaviour of the wave functions. A problem, however, exists when attempting to use 
such wave functions in descriptions of scattering from hydrogen: a credible description 
of the structure of the core is necessary to account for the full density, which is required 
to analyse nucleon scattering data. That was found in analyses of scattering from ^Li 
and ^^Li [6]. More recently, such models incorporating multiple scattering expansions 
of the nucleon-nucleon (A^A^) scattering amplitudes [7, 8] have sought to describe proton 



Modern optical potentials and the role of nuclear structure 



February 9, 2008 



1 



scattering from halos using few-body models. However, predictions from those models 
are still prone to significant changes arising from problems in specifying the density 
of the core [7]. At intermediate energies, medium modifications are essential in the 
specification of the NA optical potential (OMP) and the multiple scattering expansions, 
taken to second order [7], are only a gross approximation. A full ^-matrix specification 
of the optical potential is needed to account for those corrections [9]. 

As microscopic models of the nucleon-nucleus (NA) optical potentials are based on 
effective A^A^ interactions, these problems are overcome when models of structure which 
admit nucleon degrees of freedom are used. Herein, a description of the Melbourne op- 
tical potential [9] which accounts for scattering from both stable and exotic nuclei self- 
consistently is presented. That model has been used successfully with the shell model 
[9], Skyrme-Hartree-Fock (SHF) models [10], and the RPA [11,12]. After giving a brief 
description of the model, various results are given with a view to future experiments 
which may be done at radioactive beam facilities. Concluding remarks follow. 



THE G-FOLDING OPTICAL IVIODEL - IVIELBOURNE FORCE 

The Melbourne force is discussed in detail in a recent review article by Amos et al. [9], 
to which the reader is referred. A brief overview of the model is presented herein with 
emphasis on how structure enters into the description of elastic and inelastic scattering. 

The optical potential for NA scattering is associated with the elastic scattering chan- 
nel. Following the Feshbach formalism [13], we split the Hilbert space into the elastic 
scattering channel (P space) and non-elastic channels (Q space). The Schrodinger equa- 
tion then becomes, with P and Q projectors onto the respective spaces 



{E-Hpp) 
(E-Hqq) 



PQ 



QP 



(1) 



where Hxy — XHY. Recoupling, and taking the one-body approximation gives the 
Schrodinger equation for the projectile wave function, viz. 



[e-Ho- \V\ - (4>,, |yG^+gV| 4>,,) } \x^) = 
where g'qI = [E —Hqq -\- /e]~\ from which the optical potential is defined as 



(2) 



U = {%s\V\%s) + {^. 



■gs 



(3) 



Specification of the optical potential is a many-body problem with explicit depen- 
dence on the target ground state wave function. It is complex, nonlocal and energy de- 
pendent, through V and also Gqq. The second term is the Dynamic Polarising Potential 
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(DPP) and defines how coupling to nonelastic cliannels varies with energy. Specifically, 
such coupling may be cast into three energy regimes: 

Low energy Fori? < 10 MeV, explicit coupling to specified, discrete low-lying excited 
states of the target is necessary. A recent development in the construction of such 
models is the Multi-Channel Algebraic Scattering theory [14, 15] which is the 
subject of another presentation at this meeting [16]; 

Giant resonances Between 10 and 25 MeV coupling to the giant resonances becomes 
important [17]. One important exception are the He isotopes, for which there are 
no giant resonances; 

Intermediate and High energies At higher energies and as the level density becomes 
high, coupling to excited states may be handled implicitly by using folding models 
based on the A^A^^ g matrices for infinite matter. 

An example of the last is the Melbourne ^-folding model [9]. The model takes as its 
basis an effective A^A^ interaction obtained from the g matrices of the bare A^A^ interac- 
tion. For an incident nucleon, momentum po in collision with a nucleon embedded in 
infinite matter with momentum pi, those g matrices are solutions of the Bruckner-Bethe- 
Goldstone equation in momentum space, viz. 

where k = (po — pi ) is the relative momentum and K is the centre-of-momentum of the 
two particles. Primes denote the equivalent set of momenta after scattering. 2 is a Pauli- 
blocking operator and the energies E in the propagator contain auxiliary potentials which 
model the effects of the nuclear medium [18]. As Q and energies E are dependent on k', 
in practice they are replaced by their angle- averaged values. This has been shown to be 
a good approximation for nuclear densities above ~ 15% [19, 20], and is an important 
consideration for scattering from exotic nuclei where scattering is observed as from the 
core in the case of halo nuclei [21]. 

Once the g matrices in infinite matter are obtained, they are mapped to those for finite 
nuclei in coordinate space [9] by folding in the specified (model) density of the target. 
The mapping to coordinate space is achieved by means of a double Bessel transform 
and allows for the explicit specification of central, tensor, and two-body spin-orbit 
terms as sums of Yukawa functions. (This is also a practical consideration: the DWBA 
suite of programs [22] which are used to calculate observables require a coordinate 
space representation of the potential.) Once these effective g matrices (geff) have been 
obtained, the nonlocal, complex, OMP for scattering is defined as 

U{r,r';E) = 8{r -r')Y,ni f (p*{s)gD{r,s;E)(pi{s)ds 

i 

+ Y,ni(Pi{r')gE{r,r';E)(pi 

i 

= UD{r,E)5{r-r') + UE{r,r';E) (5) 

where D and E denote the direct and exchange terms of the effective interaction respec- 
tively. Nuclear structure information enters via the occupation numbers for each orbit 
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/. The direct term is the well-known gp form of the optical potential and is local. The 
nonlocality arises from the explicit exchange terms; neglecting such terms can lead to 
serious problems [23]. A credible model of structure is necessary in the specification of 
the optical potential. 

The single particle wave functions entering in Eq. (5) are usually assumed to be 
of harmonic oscillator (HO) form. For most nuclei this is a reasonable assumption 
and is consistent with the underlying shell model. However, for halo nuclei it is more 
appropriate to use Woods-Saxon (WS) wave functions [21] with binding energies of the 
orbits occupied by the halo set to the separation energy of the single nucleon in the halo. 

Inelastic scattering may be calculated in a distorted- wave approximation (DWA) with 
the geff as the operators effecting the transition. The transition amplitude may be 
written, with '0' and '1' denoting the projectile and bound state nucleon, respectively, 
as 

7}X'"''(^) = {xlr^ (0) I (^W l^^eff(0, l)=^oi { I ^7.M,.) (0)) } . (6) 

In Eq. (6), the distorted wave function for the projectile is denoted by and j^/qi is the 
antisymmetrization operator for the projectile and bound state nucleon. For a spin-zero 
target, one obtains after expanding the many-body wave function 







a\m\a2m.2 sj^f'r^ \ 
X (;^ir^(0)|(<Pa2(l)|Ageff(0,l)=^0l{|<Pai(l))|;^:i'"^(0))} , (7) 

where a = {l,m,j}. 



NUCLEAR STRUCTURE FACETS 

The optical model potential in the ^-folding model is a one-body operator with respect 
to the target (bound) nucleons, and so one requires specification of the one-body density 
matrix elements (OBDME), viz. 



'>aiaaJ 



= J 



'f 



Ji 



(8) 



Various models have been utilised but, for the most part, the shell model has been used 
to specify the OBDME. Others include the Skyrme-Hartree-Fock (SHE) and the RPA, 
results from both of which are presented below. Note that in specifying the OMP one 
must keep to the level of the density matrix elements as that preserves the nonlocality. 
Use of the density itself requires gross approximations to be made in the handling of the 
nonlocal exchange terms. That may be problematic [23]. 
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FIGURE 1. Differential cross sections (left) and analysing powers (right) for the elastic scattering of 
65 MeV protons from various nuclei to mass 64. 



The BonnB A^A^ interaction was used to obtain the ge/f for all results presented herein. 
All results were obtain using DWBA98 [22] from single-shot calculations; there was no 
fitting to data. The review [9] present most results obtained up to that time, and includes 
a discussion on the connection between electron and proton scattering. A subset of those 
results are discussed below, as well as those obtained since the review. 

Systematic analyses of elastic scattering across the mass range, and for several ener- 
gies, have been reported (see, eg., [9]). Results of analyses of elastic scattering differ- 
ential cross sections and analysing powers for the scattering of 65 MeV protons from 
nuclei up to mass-64 are presented in Fig. 1. Clearly the differential cross sections and 
analysing powers at 65 MeV are well reproduced by the model. Of particular note is the 
excellent reproduction of the observables' dependence with momentum transfer as one 
increases the mass. 

Fig. 2 displays the differential cross section for the elastic scattering of ^He ions from 
hydrogen at 41 A MeV (data of Lagoyannis et al. [24]) as well as the inelastic scattering 
cross section to the 2+ state. The use of WS functions, depicted by the solid line, to 
specify the density of ^He as consistent with a neutron halo gives better agreement with 
the elastic scattering data. (The results using HO functions are given by the dashed lines.) 
At these energies, the proton does not observe the halo directly but rather its effect in 
depleting the neutron density in the core to the halo. The effect is observed as a decrease 
in the differential cross section at large angles. As E2 transitions are surface-peaked, the 
halo is better illustrated in the inelastic scattering to the 2+ state as an enhancement in 
the cross section around 30°. 

The SHF model has been used to describe systematic behaviour in exotic nuclei as 
one approaches the drip lines. One may evaluate the wave functions obtained therefrom 
in analyses of reaction cross section data. We compare the results of calculations made 
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at 65 and 200 MeV using SHF wave functions with the SkX* force [25] in Fig. 3 with 
estimates by Carlson [26]. The level of agreement is quite good and the calculations 
exhibit the Coulomb shift, which is more pronounced at 65 MeV. Of note is that as one 
increases neutron number one moves away from the line of minimal isospin. 
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FIGURE 2. Differential cross FIGURE 3. Reaction cross sections 

sections for the elastic and inelastic for S, Ar, and Ca isotopes at 65 and 

scattering of ^He from hydrogen. 200 MeV. 

This illustrates a possible signature for exotic structures. 

Fig. 4 shows the proton and neutron densities for the Sn isotopes from ^'^'^Sn to 





FIGURE 4. Neutron (left) and proton (right) densities for the Sn isotopes from the SHF model. 

^^^Sn as obtained from an SHF model using the SLy4 force [27]. Those densities exhibit 
various structures as one increases mass, notably an indication of neutron halos around 
A = 150. Fig. 5 shows the differential cross sections and polarisations for a sample of 
Sn isotopes as compared to data at 40 MeV. The level of agreement between the results 
and data is reasonable and extends naturally to the variation with mass. 

Fig. 6 displays the results for elastic and inelastic proton scattering from ^osp^^ ^2], 
at 121 and 135 MeV, respectively. The RPA model was used to obtain the ground state 
and transition densities as required. The elastic scattering cross section data [28] are also 
compared to the results of an SHF calculation using the SkM* force [10]. The agreement 
with data illustrates that both the SHF and RPA models for ^ospj, fQj. ^j^g ground state are 
reasonable. Of the two models, only the RPA may be able to specify transitions densities 
and those have been used to obtain the results of the inelastic scattering. Data [29] and 
results for the transitions to the 2+ and 3^ states are compared in the right panel of 
Fig. 6. Use of the RPA allows for a self-consistent analysis of data leading to the even 
and odd parity states in ^^^Pb. That both sets of results are out of phase with respect to 
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FIGURE 5. Differential cross sections and polarisations for 40 MeV proton scattering from a sample 
of Sn isotopes. 
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FIGURE 6. Differential cross section for 121 MeV elastic (left) and 135 MeV inelastic (right) proton 
scattering from ^"^Pb. 



each other is consistent with the phase rule of Blair [30]. 

CONCLUSIONS 

We have presented a predictive model of nucleon-nucleus scattering at intermediate 
energies for which nuclear structure plays a central part. As such it is most appropriate 
for eliciting structure information of exotic nuclei from analyses of scattering of beams 
from hydrogen. Such analyses are complementary to electron scattering which, given the 
absence of electron scattering data, makes scattering from hydrogen the best (current) 
means of understanding exotic nuclei at a microscopic level. 

The model has been tested by analyses of proton scattering data across the mass 
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range and for various energies. For light nuclei, use of the shell model with appropriate 
choices of single particle wave functions allows for analyses of data out to the drip lines. 
Extension of such analyses to heavy nuclei is achieved by use of the Skyrme-Hartree- 
Fock and RPA models. In all cases, the underlying proton scattering model is consistent; 
analyses is predicated on a credible choice of structure model for the target. 

It is hoped that facilities utilising exotic beams continue scattering experiments with 
hydrogen as the target, to allow for a deeper understanding of the structures of nuclei 
out to the drip lines. Such experiments are not restricted to intermediate energies, for 
which the Melbourne force is appropriate. With the advent of a Multi-Channel Algebraic 
Scattering theory (see [16], this conference) analyses of data from low-energy facilities 
may be done with the same predictive power. 
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